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We consider a fully asymmetric one-dimensional diffusion model with mass-conserving coalescence. 
Particles of unit mass enter at one edge of the chain and coalesce while performing a biased random 
walk towards the other edge where they exit. The conserved particle mass acts as a passive scalar 
in the reaction process A -^r A —> A, and allows an exact mapping to a restricted ballistic surface 
deposition model for which exact results exist. In particular, the mass-mass correlation function is 
exactly known. These results complement earlier exact results for the A -\- A ^ A process without 
mass. We introduce a comprehensive scaling theory for this process. The exact analytical and 
numerical results confirm its validity. 



Diffusion-limited chemical processes are at the focus 
of recent research. These are dynamic systems where the 
chemical reaction time scales are short compared to those 
controlling spatial fluctuations in concentration. The lat- 
ter dominate the kinetics, in particular in low dimen- 
sions. Such processes display dynamic scale invariance 
with scaling properties that are robust, tend to be uni- 
versal, and not sensitive to many details of the actual 
dynamics at the microscopic level. Simplified models are 
therefore able to catch the essence of the process. More- 
over, some of these models are accessible to exact solu- 
tions in one dimension. 

An example of these is the one-species coalescence pro- 
cess, A -\- A ^ A Exact results for this process 
were obtained recently using the so-called inter-particle 
distribution function (IPDF) method, also known as the 
method of empty intervals ||^ . This includes several ver- 
sions of the model, with and without external input of 
particles and in the presence or absence of a diffusion 
bias along the chain [^|J|]. In this letter we address the 
fully asymmetric diffusion case where particles enter at 
one edge of the chain, and coalesce when they meet each 
other, while performing a driven random walk towards 
the other edge (I). We enhance this model by assigning 
a mass to each particle, which is preserved during each 
merging event. 

Consider a linear chain with L sites. Particles of unit 
mass enter at the left boundary x = 1, diffuse to the 
right and coalesce when they meet, and ultimately exit 
at the right boundary x = L. The diffusion of the par- 
ticles along the chain is totally biased. Choose a site x 
at random. If occupied, the particle at this site moves to 
the next site x — s- x -t- 1. If site a; -|- 1 is already occupied, 
the two particles merge, 

mx{t + 1)^0 and mx+i{t + I) = rrix+iit) + rrixit) 

(1) 



with TO the total mass of each particle. Total mass is 
conserved during coalescence. Site x = 1 is the in- 
put boundary. If chosen as an update site, it is im- 
mediately refilled by the reservoir mi{t -I- 1) = 1 and 
m2{t + 1) = m2{t) -\- mi{t) At the opposite edge, x = L, 
the particles simply fall off the chain, mi(i+ 1) = 0, back 
into the reservoir. 

The mass rides like a passive scalar on top of the par- 
ticles. It has no effect on the transition probabilities. We 
could describe the process just as well in terms of only oc- 
cupation numbers, Cx = 0, 1, for empty or occupied site. 
The latter is the fully asymmetric A -\- A ^ A process. 

The mass is a useful parameter. It allows an exact 
mapping onto an exactly soluble restricted ballistic sur- 
face deposition (RED) model Consider a one di- 
mensional interface, as shown in Fig .|l|. The steps m 
the downward staircase can take any magnitude = 
0, 1, 2, • • •. During each time step, one of the columns at 
x -|- i is chosen at random, nix particles are deposited 
onto it such that the entire step fills up and the step at 
X -\- 1 grows to nix+i + mx. Exact results for this RBD 
model have been obtained earlier using a generating func- 
tion approach . We can reinterpret those exact results 
in framework of the asymmetric A + A A reaction 
process. 

Starting from the master equation, closed form recur- 
sive equations of motion are obtained for the mass dis- 
tribution along the chain M{x,t) = {mx{t)) and also (as 
discussed later) the two-point mass correlations. M{x,t) 
obeys the relation 
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M{x, t+l) = {l- -)M{x, t) -f -M{x - 1, t) (2) 

ij Lj 

with boundary condition A/(l,t) = 1 (see eqs.(8)-(15) in 
reference Q for details). This implies that in the station- 
ary state the mass distribution is uniform, M{x) — 1, 
and there is a direct link between particle concentra- 
tion C(x) = (cx) and the average mass carried by each 
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particle, M{x). They are exactly related as M{x) = 
M{x)/C{x) where M{x) measures the average mass over 
occupied sites only. Particle coalescence creates increas- 
ingly heavier particles in the chain towards the right. At 
the same time they become more sparsely spaced, be- 
cause on average the mass remains constant along the 
chain. 

In the alternate RBD surface representation, this 
means that although the step heights increase along the 
surface from left to right, the steps occur at correspond- 
ingly larger intervals, such that the stationary state av- 
erage slope of the staircase remains constant. 

This dynamic process has a peculiar hierarchical struc- 
ture. Consider the path of one unit of mass. It performs 
a biased random walk, but completely uncorrelated from 
all other masses. Similarly, two specific mass units per- 
form completely independent random walks, until the 
moment they meet. From there on they move randomly 
but as a bound pair. Again, all other particles play no 
role. 

This hierarchical property suggests the following scal- 
ing theory. Consider a specific unit of mass. It moves to 
the right with average velocity v = 1. The standard devi- 
ation in its position is proportional to the square root of 
the time and distance traveled, i.e. Ar ~ a;^ . While dif- 
fusing to the right this particle merges with other masses. 
The total amount of mass it sweeps up is expected then 
to be of order Ar, i.e., that the average mass of occu- 
pied sites scales as M{x) ^ x^ . The amount of swept-up 
mass does not depend on the particle concentration C{x), 
because although the merging events reduce in frequent, 
they increase in size. However, this increasing lumpiness 
will show in increasing statistical fluctuations along the 
chain in Monte Carlo simulations. 

Next, we predict that the above random walk expo- 
nents are exact. This presumes the absence of intri- 
cate correlations between the particles. The hierarchical 
structure of the equations of motion, and the diffusion 
equation structure of recursion relations in both exact 
solution methods suggest this prediction. In the follow- 
ing we demonstrate the validity of this scaling theory 
from both numerical and analytical exact results. 

Hinrichsen et al. Q obtained the exact stationary state 
particle concentration. It scales as 
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This agrees with our scaling theory, because of the exact 
relation M{x) = M{x)/C{x). The scaling argument pre- 
dicts that M{x) grows as x^ , while eq.(|^) implies that 
M{x) ^ 1. 

Figure ^ illustrates the above numerically. It shows 
the time evolution of the particle density at various val- 
ues of X starting from the uniform initial state where 
every site is occupied by one unit of mass particle. We 



averaged over 10 independent Monte Carlo runs. Ini- 
tially the curves coincide, until the time when particles 
from the input edge reach that specific site. As expected 
this crossover time scales linearly with a:, i.e., with the 
uniform average particle velocity along the chain. The 
slope of the initial curve is C{x,t) ^ t^2^ and the sta- 
tionary state plateau values scale as C{x) ~ a;^ 2 ^ both in 
accordance with the scaling theory. Notice also that the 
statistical fluctuations increase with the distance from 
the source x. 

The mass auto-correlation function, = (m^) — 
(toj;)^ measures these fluctuations, and is a special case 
of the two-point correlator discussed below for which we 
have an exact solution from the mapping to the RBD 
model. At large x the fluctuations grow as Wmix) ~ y^, 
in accordance with our random-walk based scaling the- 
ory. We checked numerically the fluctuations in the par- 
ticle density. Fig.^ shows that they decay as Wp ~ A/y^, 
again consistent with the scaling theory. 

Consider the particle-particle correlation function, 
gp{x,r) = {cxCx+r) ~ {cx) {cx+r) and the mass-mass corre- 
lation function, gm{x,r) = {m^m^^r) — (itt-x) {"m-x+r) ■ As 
far as we know, no exact results are available for gp{x, r), 
from the method of empty intervals, although it seems 
within reach. On the other hand, the two-point mass 
correlator Gm{x,y,t) = {mx{t)my{t)) was shown in j|] 
to obey in the stationary steady state the recursion rela- 
tions 

Gm{^,x) = G'„i(l, a; — 1), 

Gmix,x) = 2Gmix - l,x) + Gmix - l,X - 1), 

Grn{x — 1, a:)= — G„i(a; — 2, a;), 

G,n{x,y) = i[G,„(a;- 1,?;) + G„i(a;,y- 1)], (4) 

where |a; — ?/| > 2. This set of coupled equations can be 
solved exactly, 

Gm(li x) = 1, 

G,„(a;,a;) - 1 + 4(a; - l)Z2(,_i) (0), 
Gm{x - l,a;)= Z2{x-2){^), 

r-l v-2 

G„i{x,y) ^^Z2{y-2){n) = l-^Z2(y-2)in), (5) 

n—O n—r 

with r = (y — a;) > 0, and 
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{2p-n)\ ,1.2p^„ 



p\{p-n)V2' 



(6) 



which is related to the probability that a random walker 
returns to its starting point exactly n times up to 2p 
steps, see for details. 

We expect the following scaling form for both gm and 
gp in the stationary state: 
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5o(a:,r) =6-2-=g„(6-ix,6-i/2r), (7) 

with b an arbitrary scale factor. The distance r from the 
source plays the role of time in our diffusion-coalescence 
type scaling argument. Therefore r should scale with re- 
spect to the correlator distance x a.s r ^ y/x. The other 
exponents, Xp and x^, follow by power-counting. They 
must be the dimensions of the particle and mass con- 
centrations, C{x) and M{x), which are equal to Xp = i 
and Xm = according to the previous discussion. We 
perform numerical simulations for gp, and exact enumer- 
ations of the exact formula for (/,„. Fig.^ and Fig.||show 
the scaling functions fp and /„, defined as 

gm^ f,n{r/Vx), (8) 

The data collapses perfectly and thus confirms the valid- 
ity of the scaling relations, eq.(0). 

Both scaling functions vanish at large R = rj^/x, as 
they should. At small R they are both linear in the scal- 
ing variable. It is easy to evaluate the exact recursion 
relations for Gm in the limit of large x and fixed small r. 
This yields ~ -l-f -^R. 

The particle correlator scaling function is also linear at 
small R. We find numerically that 

fp{R)^-a + bR (9) 

with a = 0.318(1) and b = 0.268(6) and we suspect that 
the overall factor a = I/tt. 

The above numerical and exact results are all in full 
agreement with the simple scaling picture where we treat 
the particles as performing free independent biased ran- 
dom walks before they merge. None of the scaling expo- 
nents differ from their diffusion values. The final verifi- 
cation for the validity of this intuitive explanation is the 
shape of the stationary state mass distribution function, 
plrrijx), i.e., the probability to find a particle of mass m 
at a site x 3> 1. Our scaling theory presumes that this 
distribution behaves in the same manner as the proba- 
bility for a specific tagged particle to reach site x with 
mass TO. Each tagged particle follows an independent bi- 
ased random walk and its mass grows proportional to the 
spatial fluctuations about its average path. 

Fig.^ shows numerical results for p(to, x) at system 
size L=512 for various values of a;=16, 64, 256, and 
512. It is obtained from 10 independent runs, each 10^ 
Monte Carlo steps long. The data is plotted in terms of 
y/xp{m,x) versus the scaling variable ^ = m/^/nx and 
collapses well onto Gaussian with a linear prefactor, 

p{m,x)^A — exp(-i3TO^/a;). (10) 
x 

The parameters A and B are predetermined by the nor- 
malization condition and the scaling of the particle den- 
sity (see eq.(|)) 
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^ p{m,x) = 1, 

m=0 

oo 

M{x)^'^mp{m,x) = l/C{x), (11) 

m=0 

which yields A — ^ and B ~ j. The drawn curve in 
Fig.^ corresponds to this Gaussian with the above values 
of A and B. Eq.(|l^) is the simplest Gaussian form con- 
sistent with the requirement that p{m, x) vanishes in the 
limit TO — !■ 0. According to our intuitive picture, p{m, x) 
is related to the probability that a biased random walker 
with drift velocity = 5, makes an excursion of size to 
from its average path before reaching site x, irrespective 
of (i.e., averaged over all) the starting times at site a; = 1. 

In conclusion, we present new exact results for the 
A + A ^ A type coalescense process by introducing mass 
to the particles, and exploring the exact mapping to a 
surface deposition model, the so-called RBD type surface 
growth model. In addition, we propose a scaling theory 
based on the assumption we can treat the particles as 
performing free independent biased random walks before 
they merge. All scaling exponents should then take naive 
diffusion values. The above numerical and exact results 
are all in full agreement with this random walk type scal- 
ing. It appears therefore that the scaling properties of 
A + A ~^ A type dynamics are now fully understood, 
and actually, in the final analysis, are predictable from 
random walk considerations only. 
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FIG. 3. fluctuations in the particle concentration from the 
auto-correlation functions Wp{x) = {c?{x)) — {c{x))'^ at chain 
length L = 512. The slope of the drawn line is -1/2. 
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FIG. 1. Surface representation of the A + A A process. 
The step heights along the stairway represent the mass of 
the particles. Each deposition process fiUs-up the entire step, 
such that the adjacent step has a new height equal to the sum 
of the two. 




R=[r/x^'^] 

FIG. 4. The scaling function for the particle-particle cor- 
relator gp, from 10 independent Monte Carlo runs each aver- 
aged over lO'' time steps at chain length L — 512 for various 
distances from the source, r — 1(+), 2(x), 4(o), 8(A), 16(v)j 
and 32(o). The drawn line is obtained from eq.(p|). 
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FIG. 2. The time dependent particle concentration {c{x)) 
as function of time at various values of the distance from the 
source, x = 16(x), 32(o), 64(.), 128(A), 256(v), and 512(o) 
at chain length L = 512. The slope of the drawn line is -1/2. 
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R=[r/x^'^] 

FIG. 5. The scaling function of the mass-mass correlator 
Qm, from enumeration of the exact recurrence relations, at 
chain length L = 512 for various distances from the source, 
r = 2(x), 4(o), 8(A), 16(v), and 32(o). The drawn line is the 
tangent at small R, see between eq.M and eq.(0). 
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FIG. 6. Numerical Monte Carlo scaling plot of the station- 
ary state distribution p{m, x), the probability to find a parti- 
cle of mass m at site x, for various x = 16(x ), 64(o), 256(«), 
and 512(A). We plot ^/xp{x,m) versus ^ = m/v^ra, the 
Gaussian form of eg. ([lo|) . 
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